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Abstract 

We report a quasi-exact power law behavior for Ising critical tem- 
peratures on hypercubes. It reads J/fc^Tc = i^o[(l — ^/d){q — 1)]" 
where Kq = 0.6260356, a = 0.8633747, d is the space dimension, q 
the coordination number [q = 2d), J the coupling constant, kg the 
Boltzman constant and Tc the critical temperature. Absolute errors 
from available exact estimates (d = 2 up to d = 7) are always less 
than 0.0005. Extension to other lattices is discussed. 
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In ref. [1] we showed that both site and bond percolation thresholds obey 
a universal power law, 

p, = Po[{d-l){q-l)]-^d' (1) 

where d is the space dimension, q the coordination number, po ^-^d a are 
constants. For site dilution 6 = while 6 = a for bond dilution. 

From a Ln — Ln plot, all known percolation thresholds are found to align 
on two straight hnes each one thus defining a universal class characterized 
by a set of parameters {po] a}. One class includes two-dimensional triangle, 
square and honeycomb lattices, with {po = 0.8889; a = 0.3601} for site dilu- 
tion and by {po = 0.6558; a = 0.6897} for bond dilution. Two-dimensional 
Kagome and all lattices of cubic symmetry (for d> 3) constitute the second 
class with {po = 1.2868; a = 0.6160} and {po = 0.7541; a = 0.9346} for sites 
and bonds respectively. At high dimensions a third class for hypercubes (sc 
and fee) is introduced to recover the infinite Caylcy tree limit [1]. We noticed 
that some Ising critical temperatures were found to obey Eq. (1) within the 
bond scheme {b = a). 

Very recently Hackl and Morgenstern suggest a connection between on the 
one hand bond percolation thresholds and, on the other one, a combination of 
Ising critical temperatures and energies [2]. Agreement with available data 
are within few percents. 

In this note we report a quasi-exact power law formula for Ising critical 
temperatures on hypercubes, 

K^^Ko[{l-^){q-ir\ (2) 

where K^. = J is the coupling constant, ks the Boltzman constant, 

Tc the critical temperature, q = 2d, Kq = 0.6260356 and a = 0.8633747. 
Absolute errors are always less than 0.0005. Apphcation of Eq. (2) to other 
lattices is discussed. 

Fig. (1) shows a Ln — Ln plot of respectively K^^ and [(1 — l/d){q — 1)] 
for all Ising hypercubes from d = 2 up to d = 7. Data arc exact estimates 
taken from refs. [3, 4]. All points (small squares in the Figure) align on a 
single straight line which demonstrates Ising critical temperatures (Kc) obey 
a power law of the form of Eq. (2). 

Characteristics of that straight line are determined using exact estimates 
from d = 2 up to d = 6 not including d = 7. We found Kq = 0.6260356 
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and a = 0.8633747. Associated values of Kc are reported in Table (1). The 
agreement with exact estimates is indeed remarkable. Defining A = Kc — 
where are exact estimates [3, 4] and Kc are from Eq. (2), absolute errors 
|A| are found to be always less than 0.0005. 

Including d = 1 yields Kq = 0.6269574 and a = 0.8647815 which still 
is excellent but a little bit less accurate. The error A is then reduced for 
d = 7 down to A = +0.00024 but goes up to A = +0.00084 for d = 2 (see 
Table 1). This sensitivity to the case d — 7 sheds light on our conjecture of 
a crossover to a new universality class for percolation thresholds above d — 6 
[1]. Moreover Eq. (2) with a ^ 1 does not have the d oo asymptotic 
Caylcy tree limit Kc = 1/q hinting to some crossover at high dimension. 

In Fig. (2) Ising critical temperature exact estimates for non-hyper cubic 
lattices honeycomb, kagome, triangular at d — 2, and diamond, bcc and fee 
at d = 3 are included. Data arc seen not to sit exactly on the hypercubic 
straight line. Discrepancies are between A = —0.031 and A = 0.0082 (see 
Table 2). It hints some correction to q = 2d should be accounted for in Eq. 
(2). We found for instance that rescaling the variable (1 — l/d){q — 1) by 
{q/2d)^^^^ reduces discrepancies for non-hypercubic systems. However this 
point needs more investigation. 

To complete our presentation we have also determined Kq = 0.6247099 
and a = 0.8606241 using only d = 2 and d = 3 data. Associated errors A* 
are shown in last column in Tables (1, 2). 

To conclude we have found a quasi-exact formula (Eq. 2) to yield Ising 
critical temperatures on hypercubic lattices. At this stage it worthwhile to 
stress that a posteriori our results provide ground to technical efforts made 
to determine critical temperature numerical estimates up to more than six 
digits [3, 4]. Otherwise we would have conclude wrongly Eq. (2) is exact as 
seen from associated errors in Table (1). 

More estimates at d — S and d — 9 are now needed to elucidate the high 
dimension crossover to the d ^ oo asymptotic Caylcy tree limit. Last but 
not least, it is interesting to note [(1 — exp(— 2i^rc)] does not obey a power 
law in the variable [(1 — l/d){q — 1)]. 
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Figure captions 

Figure 1: Ln — Ln plot of versus for hypercubes from d — 2 

up to d = 7. Small squares are exact estimates from [3, 4]. The straight line 
is from Eq. (2). 

Figure 2: Ln — Ln plot of K^^ versus for hypercubes and other 

lattices from d = 2 up to d = 7. small squares are exact estimates from [3, 
4]. The straight hne is from Eq. (2). 
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Dimension 


Lattiee 


Q 




Kc 


A 


A* 


d = 2 


Square 


4 


0.4406868 


0.44112972 


+0.00044 







se 


6 


0.2216544 


0.22139040 


-0.00026 





d = 4: 


sc 


8 


0.14970 


0.14956561 


-0.00013 


+0.00023 


d = 5 


sc 


10 


0.1139291 


0.11386752 


-0.00005 


+0.00032 


d = 6 


sc 


12 


0.092295 


0.09228595 


+0.00014 


+0.00051 


d=7 


se* 


14 


0.077706 


0.07809958 


+0.00039 


+0.00075 



Table 1: Ising critical temperatures from this work Kc compared to "exact 
estimates" taken from [3, 4]. A = A'c — K^. * means not included to 
determine Ko = 0.6260356 and a = 0.8633747. A* corresponds to Kq = 
0.6247099 and a = 0.8606241 using only d = 2 and d = 3 data. 



Dimension 


Lattice 


q 


K! 


Kc 


A 


A* 




Honeycomb* 


3 


0.6584788 


0.62603564 


-0.032 


-0.034 


d = 2 


Kagome* 


4 


0.4665661 


0.44112972 


-0.025 


-0.026 




Triangular* 


6 


0.2746531 


0.28381002 


+0.009 


+0.009 


d = 3 


Diamond* 


4 


0.36982 


0.34411006 


-0.026 


-0.026 




bcc* 


8 


0.15740 


0.16557529 


+0.008 


+0.008 




fee* 


12 


0.10209 


0.11207783 


+0.010 


+0.010 



Table 2: Ising critical temperatures from this work Kc compared to "exact 
estimates" taken from [3, 4]. A = Kc — K^. * means not included to 
determine Ko = 0.6260356 and a = 0.8633747. A* corresponds to Kq = 
0.6247099 and a = 0.8606241 using only hypercubic d = 2 and d = 3 data. 
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